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1. INTRODUCTION 

In power systems, the hybrid active power filter (HAPF) is considered as the most effective solution 
in harmonic filtering and reactive power compensation [1-3]. However, the working efficiency of HAPF 
depends on many factors such as: correct determination of the parameters [4], control method [5-8], harmonic 
current detection method [9, 10], control strategy [11, 12] and stabilization of the DC-bus for the inverter [13, 
14]. In which, the correct calculation of the parameters of the HAPF is the most important factor, it decided to 
work efficiency of HAPF. 

The structure of a HAPF includes: power circuit part and control circuit part. Therefore, the 
parameters that need to be designed for HAPF are the parameters of the above circuit parts. Until now, the 
multi-objective optimization studies of HAPF can be summarized as follows: The application of the genetic 
algorithm to multi-objective optimization design for passive power filter (PPF) is proposed by [15]. Another 
algorithm is also commonly used for HAPF design, called the Particle Swarm Optimization (PSO) algorithm. 
The PSO algorithm is used to design an APF in a four-wire three-phase system in case of balanced and 
unbalanced loads but only consider to design for PPFs and optimization for the APF but not for the HAPF [16- 
19]. Another study to multi-objective design on HAPF is using an ant colony and bat algorithm [20, 21] but 
only aslo multi-objective optimization design for PPFs. Ahmed Faheem Zobaa has used the Fortran Feasible 
Sequential Quadratic Programming algorithm to solve the multi-objective optimization problem for HAPF 
with the aim of finding PPF parameters [22]. 

In summary, the previous research on multi-objective optimization for HAPF was the only multi- 
objective optimization design of power circuit parts (parameters of the PPF). Meanwhile, the optimized design 
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for the parameters of the control circuit part has not been studied and the calculation of the above parameters 
without considering the condition of the system stability. 

To overcome this disadvantage, in this paper, the Routh’s stability standard is used to check the stable 
condition of the HAPF system. The Vortex search is a new optimization algorithm, that is effecty in solving 
global optimization [23, 24]. Then a new multi-objective optimization design, based on the Vortex Search (VS) 
algorithm is proposed to determine the best set of parameters for the HAPF. The proposed algorithm has the 
advantage of finding a quick solution with little loops. The results achieved will be global optimization such 
as the minimum of the total harmonic distortion (THD) of the supply current and source voltage and satisfy the 
stable conditions. 

The structure of the paper is divided into five parts. Part 1 is an overview introduction of the issues 
that need to be investigated. Part 2 is a new multi-objective optimization design method for Hybrid Active 
Power Filter. The simulation results and discussion are presented in Part 3, and the conclusions are summarized 
in Part 4. 


2. NEW MULTI-OBJECTIVE OPTIMIZATION DESIGN METHOD FOR HYBRID ACTIVE 
POWER FILTER 
The single-phase equivalent circuit of the selected Hybrid Active Power Filter topology is shown in 
Figure 1. It consists of five parts: U; and Z; are the source voltage and the resistance of the source. C; and L; 
are capacitance and inductance at the fundamental frequency. Cy is the added capacitance to compensate for 
reactive power. Lo is the output filter of the voltage source inverter. 
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Figure |. Single-phase equivalent circuit of the selected hybrid active power filter 





Starting from the load current iz (ira, iz, inc), We USe p-q instantaneous power theory [9], [10] to 
separate the harmonic component of the load current as izna, itn, itn. These components will be chosen as the 
reference components and they are compared to the compensate current components into the grid i, of HAPF 
are ica, ich, aNd icc. The error of the above comparison will be taken through the PI (Proportional - Integral) 
controller, through pulse width modulation to control ON-OFF switches of the inverter to generate compensate 
signals on the grid. Thus, the inverter will be controlled to produce a signal similar to the harmonic current 
signal generated from the load. 

The parameters of both the power circuit part and the control circuit part need to be determined in the 
model in Figure | including passive circuit Cr - C; - L; in which branch C; - L; resonates at the fundamental 
frequency, output filter Lo, voltage bus-DC of inverter Ua, parameters of PI controller are K, and Kj. Currently, 
these parameters are usually determined based on experience and local calculations, regardless of the stability 
of the system. Therefore, this paper proposes a new multi-objective optimization design method that has 
considered the stability of the system. 


2.1. Constraints and Objective Function 
System stability constraint: Based on the transfer function of the HAPF system in [1], we have the 
characteristic equation of the control transfer function as follows: 





D(s) =a,8° +4, +4,5' +a, +a,8° +a," (1) 
According to Routh’s stability standard, the following conditions must be satisfied. 
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ba, —a,b, >0 (2) 
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inv 








Constraints on resonance conditions in PPFs: The LZ and C parameters in a branch must be resonated 
at a certain frequency. 


he ees (3) 


Constraint of passive power filter values: passive power filter values must be positive and meet system 
stability and resonance conditions. 


0< (R,L,C)) & E airy Sean, S09 (4) 


the values of Rmax, Lmax, ANd Cmax are determined under stable conditions (2). 
PPFs must be compensated with a maximum capacity but not exceed the required maximum limit 


Q nin < Q. < Q sax (5) 
constraint of DC-bus voltage value: 
0 < Uae < Ude-max (6) 


Controller parameters constraints: Controller parameters must be positive and satisfy system stability 
conditions (2). 


O<Ky < Konae ¢ 0 < Kp< Kins (7) 


according to [5], the value of K, and Kjis too small or too large to cause system instability. 
Objective function: we consider the main objective as follows: 


F = min(THDi, + THDu, ) (8) 


2.2. Vortex Search Algorithm and Its Application for Designing Hybrid Active Power Filter 

The Vortex search algorithm is first proposed in [23, 24]. Let us consider a two-dimensional 
optimization problem. In a two-dimensional case, a vortex pattern can be modelled by a number of nested 
circles. Here, the largest circle of the vortex is first centered on the search space, where the initial center can 
be calculated using (9). A description of the VS algorithm is also provided in algorithm 1. 


__ upperlimit+lowerlimit 
0 
2 





(9) 


max ( upperlimit )— min(lowerlimit 
h YO, = ( is I ( ) (10) 
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The solutions that exceed the boundaries are shifted into the boundaries, as in (11). 


rand (upperlimit' —lowerlimit' ) +lowerlimit';s. < lowerlimit' 
s= 5, ;lowerlimit' < s, < upperlimit' (11) 


> upperlimit' 


rand (upperlimit' —lowerlimit' ) +lowerlimit';s\ 


The initial radius ro can be calculated with (10). Because ao =1, the resulting function value is 
1(x)-gammaincinv (x, ao) ~1, which means ro ~ dy as indicated before. 


1 
r= o.( +) gammaincintx, a) “ 


by methods of (12), the initial radius value rp can be computed as rp~10. Radius at each iteration is calculated 
by using (13). 


1 
r= o.(+) .gammaincinv(x, a,) or 
x 


In this algorithm paramater a, of the inverse incomplete gamma function defines the resolution of the 


search. By equally sampling values within [1, 0] interval at a certain step size, the resolution of the search can 
be adjusted. For this purpose, at generation t, a is obtained by using (14). 


t 
a4=a,- 
‘  "° MaxInt 





(14) 


where a, is selected as a,=1 to ensure full coverage of the search space at the first iteration, 


t is the iteration index, and MaxI/tr represents the maximum number of iterations. 
Vortex Search algorithm flowchart for designing HAPF is shown as in Figure 2. 





Algorithm 1. 
Inputs: Initial center gis calculated by using (9) 
Initial radius rpor dis caculated by using (10) 
Best fitness so far f (Spest) = inf 
Output: Best solution Spest 
t= 0; 
Repeat 
/* Generate candidate solutions by using Gaussian distribution around the center with a standard deviation 
(radius) 7,*/ 
Generate C,(s); 
If exceeded, then shift the C,(s) values into the boundaries as in (11); 
/* Select the best solution from C,(s) to replace the current centerp,*/ 
s' = select(C;(s)); 
if (8) < f (pest) 
Sbest = 5 


f (Svest) = f(s) 





else 
Keep the best solution found so far Spese; 
End 
/* Center is always shifter to the best solution found so far*/ 
Me+1 = Sbests 
/* Decrease the standard deviation (radius) for the next iteration*/ 
T141 = Decrease(r;); 
t=t+1; 
Until the maximum number of iterations is reached 
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Figure 2. Vortex search algorithm flowchart for designing HAPF 


Application of the VS Algorithm for Hybrid Active Power Filter Design is described as follows: The 
algorithm begins by entering the upper and lower limit values of the parameters to be searched based on the 
constraints in the formulas from (3) to (7), then these limits will be included in the VS algorithm to find a set 
of parameters to be searchedX (Cp, Cy, L4, Lo, Vac, Kp, Ki). Next, these parameters will be checked for stable 
conditions using formula (2), if not stable condition, return to VS algorithm to create a new set of parameters, 
if stable condition is satisfied, these parameters will be setup into the HAPF simulation model to find the 
THDi;% and THDu;% values. 

The value of THDi,% must be less than ¢, and THDu,% must be less than ¢,, according to IEEE 519 


standard [25], the value of ¢, is less than or equal to 5% and ¢, is less than or equal to 5%. Therefore, to satisfy 


€, and €, the objective function F must be less than 5%, then stop and print out the results. Vortex Search 
algorithm flowchart for designing HAPF is described as in Figure 2. 


3. SIMULATION RESULTS AND DISCUSSION 

To demonstrate the effectiveness of the proposed design method. Let's consider, let us consider a 
HAPF model as shown in Figure 1. Three-phase source voltage is 380V-50Hz. THD of the load current is 
shown in Figure 3. 

According to [1], the parameters of the HAPF system with traditional design method are calculated 
in Table 1. In addition, the controller parameters are randomly selected with the best parameters of K, = 30 and 
K; = 0.1, switching frequency f= 10kHz. 
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Figure 3. THD of the load current 


The parameters in Table 1, we obtain the waveforms of the HAPF system shown in Figure 4. From 
Figure 4, we see that: THD of the supply current i; decreases from 32.77% to 2.1%, THD of the source voltage 
us is 1.69%, while the reactive power decreased to 140var from 2359var, that is, the capacity compensated by 


PPFs is 2219var. THD of the supply current i; shown in Figure 5. THD of the source voltage u; is shown in 
Figure 6. 


Table 1. HAPF Parameters by Traditional Design Method 


Cr (UF) Ci (uF) Li(mH) Lo(mH) Uae (V)__THDI,(%) __ TH Du; (%) 
116.8 349.2 29.77 0.2 535 2.1 1.69 
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Figure 4. Waveforms of HAPF with traditional design method 
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with the proposed design method, the best parameters of HAPF after 10 iterations are shown in Table 2. From 
the parameters in Table 2, we obtain the waveforms of the HAPF system shown in Figure 7. 


Table 2. HAPF Parameters with the Proposed Multi-Objective Optimization Method 





Cr Cc, Ly Lo Uae K, K; THD i, THD u, 
(HP) (HF) (mH) (mH) (V) (%) (%) 
136.67 301.22 26.6 1.01 674.13 94.95 0.65 1.14 1.28 
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Figure 7. Waveforms of HAPF with the proposed design method 


From Figure 7, we can see that: THD of the supply current i; decreases from 32.77% to 1.14%, while 
the reactive power decreased to 90var from 2359var, so the compensation power by PPFs is 2269var. THD of 
the supply current i, in steady-state is shown in Figure 8, THD of the source voltage u; in steady-state is shown 
in Figure 9. 


Fundamental (50Hz) = 226.2 , THD= 1.14% Fundamental (50Hz) = 308.7 , THD= 1.28% 


























100 100 
® so} a 
2 2 
= =] 
ire irs 
‘ 4ol S 
fe = 
D D 
& 20; & 
0 
0 200 400 600 800 1000 0 200 400 600 800 1000 
Frequency (Hz) Frequency (Hz) 
Figure 8. THD of the supply current i, Figure 9. THD of the source voltage us; 


According to IEEE Std. 519-1992 Recommended Practices and Requirements for harmonic control 
in electrical power systems [25], then total harmonic distortion of current is less than or equal 5%. THD of the 
supply voltage is less than or equal 5%. From the simulation results, we can see that: the proposed design 
method is more effective than the traditional design method in reducing the total harmonic distortion of the 
supply current and source voltage and meet IEEE Std. 519-1992. 
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4. CONCLUSION 

The paper has proposed a new multi-objective optimization design algorithm for HAPF based on the 
Vortex Search algorithm. This method can determine the parameters of both the power circuit part and the 
control circuit part of the HAPF. Compared to the traditional design method, simulation results have proved 
that: the achieved results are multi-purpose, such as: minimize the total harmonic distortion of the supply 
current and source voltage, maximize the reactive power compensation into the grid and satisfying many 
constraint conditions. This research can be applied to design for all different types of HAPF. 
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